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ABSTRACT 

This paper is a sequel to an earlier paper devoted to multiple, multicycle 
nova evolution models (Prialnik & Kovetz (1995), first paper of the series), which 
showed that the different characteristics of nova outbursts can be reproduced by 
varying the values of three basic and independent parameters: the white dwarf 
mass-M^D, the temperature of its isothermal core-Two and the mass transfer 
rate-M. Here we show that the parameter space is constrained by several ana¬ 
lytical considerations and find its limiting surfaces. Consequently, we extend the 
grid of multicycle nova evolution models presented in Paper I to its limits, adding 
multicycle nova outburst calculations for a considerable number of new parameter 
combinations. In particular, the extended parameter space that produces nova 
eruptions includes low mass transfer rates down to 5 x 10“^^ Mq yr“^, and more 
models for low Two- Resulting characteristics of these runs are added to the 
former parameter combination results, to provide a full grid spanning the entire 
parameter space for Carbon-Oxygen white dwarfs. The full grid covers the en¬ 
tire range of observed nova characteristics, even those of peculiar objects, which 
have not been numerically reproduced until now. Most remarkably, runs for very 
low M lead to very high values for some characteristics, such as outburst am¬ 
plitude A > 20, high super-Eddington luminosities at maximum, heavy element 
abundance of the ejecta Zf.j ~ 0.63 and high ejected masses mej ~ 7 x 10“^ Mq. 
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1. INTRODUCTION 

The realization that Classical Novae (CN) are due to thermonuclear runaways (TNR) on 
the surfaces of white dwarf (WD) stars in close, mass-transferring binary systems (Starrheld, 
Sparks & Truran 1972) is now three decades behind us. Direct observations of novae (sum¬ 
marized in e.g. Payne-Gaposchkin (1957)) had already demonstrated that large ranges in 
ejecta velocity and metallicity occur in novae, and later interpretations of the data suggested 
that ejecta masses varied strongly from object to object. The first simulations by nova theo¬ 
rists, trying to make sense of these observations, involved much groping into various corners 
of parameter space. Early triumphs included an explanation of the UV emission and high 
bolometric luminosities of novae after eruption (Starrheld, Sparks & Truran 1976); fast and 
moderately slow nova models (Prialnik, Shara & Shaviv 1978); natural enrichment of ejecta 
in heavy elements via diffusion (Prialnik & Kovetz 1985) and via mixing from the underlying 
white dwarf (MacDonald 1983; Kutter & Sparks 1987); the realization that nuclear burning 
with EUV emission could occur without mass ejection (perhaps leading to SNIa or super 
soft X-ray sources) (Shara, Prialnik & Shaviv 1977) and novae recurring on timescales as 
short as a few years (Starrheld, Sparks & Shaviv 1987). 

The hrst decade of nova simulations crystallized the understanding that a space of three 
independent parameters controls the behavior of a classical nova eruption: the white dwarf 
mass, the white dwarf temperature (or luminosity) and the mass accretion rate from the 
companion (Shara, Prialnik & Shaviv 1980; Prialnik 1995). Increased computer power, and 
codes including better physics lead to an ever-increasingly sophisticated series of eruption 
simulations. Once multi-cycle simulations became possible, the arbitrariness of initial con¬ 
ditions could be eliminated, and systematic surveys throughout the three-dimensional space 
became a real possibility. These culminated in an extended grid of multicycle nova evolution 
models (Prialnik & Kovetz (1995), hereafter Paper I) which have been extensively used by 
researchers. 

While Prialnik & Kovetz covered virtually all of the {M^D) M) space then thought 
likely to be reached by novae, more recent work (Hurley & Shara 2003) has shown that cata¬ 
clysmic variables are likely to evolve considerably beyond the boundaries of their calculations. 
Remarkably, we hud (in the present study) that some of these systems (parameter combina¬ 
tions) still produce nova eruptions. Amongst these are models with extreme values of 
Ty/j:, and M which can give rise to ’’extreme novae” - objects with eruption luminosities, 
metallicities, and/or ejecta masses and velocities signihcantly larger or signihcantly smaller 
than any yet predicted by simulations. 

The existence of very unusual eruptive systems like the ’’Red Variable” in M31 (Mould 
2000) or V838 Mon (Bond et ah 2003) are already sufficient justihcation for extending nova 
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simulation grids to cover all conceivable cases. Each of these unusual variables has properties 
(particularly luminosity) reminiscent of novae, and we will examine possible connections in 
a future paper. But of equal interest is the physical insight that these new simulations can 
provide, as important input in trying to determine the very long term evolution of cataclysmic 
binaries, and relationships between the various subclasses of CVs. Most surprising of all is 
the non-linear behavior of some of the most extreme new nova models. Simple extrapolations 
from previous simulations give answers that may be wrong by large factors. This is because 
dominance of competing timescales and physical effects changes with movement through the 
three-dimensional nova parameter space. 

It is perhaps not too surprising that M lower than 10“^^ Mq yr“^ (the lowest values 
considered in Paper I) can occur, especially as a mass-accreting WD whittles down its donor 
companion to the brown dwarf mass range. Equally understandable is that WD core tem¬ 
peratures can drop to or below 10^ K, for WDs accreting at very low rates for many Gyr, 
whereas Paper I considered only hot WDs accreting at the lowest rate. Both these effects - 
low mass accretion rate and cold WDs - are extensively studied in the present work. What 
is less intuitively obvious is that CO WDs with masses well under the 0.65 Mq (the lower 
limit of the Paper I study) are also expected. Mass stripping via common envelope evolu¬ 
tion can produce a naked helium star, which evolves to the helium giant branch and then 
becomes a low-mass CO WD, as illustrated by the example of a long-term binary evolution 
track summarized in Table 1 (Hurley, private communication). The existence of low-mass 
(< 0.50 Mq) WDs resulting from close binary evolution was also discussed by Marsh et al. 
(1995) and Althaus et al. (2004) (and references therein). The behavior of such low-mass 
CO WDs accreting hydrogen-rich material is also examined, for the first time, in the present 
work. 

The methods of numerical computation, the grid of results, and the range of nova 
parameter space covered are discussed in § 2. In § 3 we discuss the existence of a conhned 
region of the three-dimensional parameter space within which conditions for nova outbursts 
are satished. A brief summary and conclusions are given in § 4. 

2. NUMERICAL COMPUTATIONS 

The hydrodynamic Lagrangian stellar evolution code by which the present study was 
performed is the code presented in Paper I. It includes OPAL opacities, an extended nuclear 
reactions network comprising of 40 heavy element isotopes and a mass-loss algorithm that 
applies a steady, optically thick supersonic wind solution (following the phase of rapid ex¬ 
pansion). In addition, diffusion is computed for all elements, accretional heating is taken into 
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account and convective fluxes are calculated according to the mixing length theory. Initial 
models were prepared for four WD-mass values and three temperatures by cooling WD mod¬ 
els from higher temperatures. Each nova model was followed through several consecutive 
outburst cycles in order to eliminate the effect of the initial configuration. One typical cycle 
was then chosen as representative of each parameter combination. 

The original grid presented in Paper I consisted of 64 parameter combinations: four 
Mwd values—0.65,1.00,1.25 & 1.40 Mq, three Ty/j:, values—10,30 & 50 x 10® K and five 
M values—10“® through 10“^® Mq yr“^, and four calculations for the highest temperature 
value with M of 10“^^ Mq yr“^. In the present work we have extended the grid to lower 
accretion rates, down to M = 5 x 10“^^ Mq yr“^ for all values of WD mass and temperature 
and we have added calculations for M\yd = 0.4 Mq, amounting to about 30 new parameter 
combinations in all. 

Tables 2 & 3 list the results of our complete grid of parameter combinations. The tables 
are presented in a manner similar to that of Paper I: Table 2 lists properties that are related to 
the accretion phase and the onset of the outburst, whereas Table 3 presents characteristics of 
the outburst itself. The properties displayed in Table 2 are: the accreted and ejected masses 
^acc, mej] the helium mass fraction of the envelope and the ejecta Y^nv, Yej] the heavy 
element mass fractions Zenv, Zej] and the maximum temperature attained in the burning 
shell Tmax- The characteristics displayed in Table 3 are: the maximal expansion velocity Vmax, 
its average over the whole mass-loss phase Vav, the maximal luminosity attained during the 
outburst L^^rnax, the amplitude of the outburst in magnitudes and three typical timescales— 
the time of decline of bolometric luminosity by 3 mag. t^^boi, the duration of the mass-loss 
phase tm-i, and the recurrence period of the outbursts Free = mace/M. 

Some parameter combinations yield peculiar results, such as extremely high (highly 
super-Eddington) luminosities, which will require further investigation; for a few runs the 
determination of timescales, such as ts^boi and tm-i may not be accurate, due to fluctuations 
of the light curve. For example, the combination 125.30.13 (standing for 1.25 Mq WD, with 
Twd = 30 X 10® K and M = 10“^^Mq yr“^) produces a sharp quick peak at outburst, 
resulting in very small t^^boi and tm-i (Table 3), followed by a second increase in luminosity 
and temperature, but one that does not eject more material. The values stated in the table 
relate only to the first quick rise. 

One of the most important characteristics of a classical nova outburst is the time of 
decline of the visual magnitude, which serves as a basis for the ’’speed class” classification 
of novae and on which the role of novae as standard candles is based. In Fig. 1 we show 
the distribution of the time of decline (which is of the order of tm-i as presented in 
Table 3) versus the outburst amplitude A for 79 observed galactic novae from the Duerbeck 
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catalog (Duerbeck 1987), together with results of the sub-set of 75 models that produced 
nova outbursts. Three points should be noted. First, the general tendency of decreasing for 
increasing A is clearly seen. Secondly, the model results provide full coverage of observations 
in the [A, t^] plane. In this respect, it should be noted that the density of points on the plane 
has no signihcance, since the calculated models represent parameter combinations (they are 
not based on a population synthesis model). Finally, we predict the existence of a class of 
low amplitude novae (with A ~ 7) with decline times ranging from 7 days to over a year. 
We suggest that these objects have either not yet been observed, or remain characterized as 
variables of some other type. 

The effect of each of the three input parameters on the resulting nova characteristic was 
already discussed in Paper I and also by Schwartzman et al. (1994). The outcome is not 
always self-evident and the trends are not necessarily monotonic. To illustrate the behavior 
of some of the more interesting characteristics, we plot in Fig. 2 the model results over the 
{Mwd, M) plane for a given T^d value. Panel 2a displays the maximum luminosity attained 
Li,max together with the LEdd surface, enabling the identihcation of super- and sub-Eddington 
luminosities. The lowest sub-Eddington luminosities are obtained for the low Mwd, high 
M corner, identihed as the region of least-powerful nova outbursts. Panel 2b displays the 
ejected mass rriej (on a logarithmic scale). The highest ejected masses are obtained for the 
lowest WD masses through the entire range of accretion rates. This is because the accreted 
masses necessary to trigger a TNR are higher for lower mass WDs (see § 3) whereas at 
the same time, the weaker gravitational potential enables more massive ejecta, despite the 
typically lower outburst intensity. In panel 2c we show the maximum expansion velocities 
Vmaxi which exhibit a similar behavior to that of the maximum luminosity surface. While 
the low My/D, high M domain is that of weakest novae, the high Mwd, low M corner can 
be identihed as the domain of the most powerful outbursts. The time of bolometric decline, 
as shown in panel 2d (the trend being very similar to that of tm-i, which is not displayed 
here), clearly exhibits, in conjunction with hgures 2a and 2c, the correlation between decline 
time (speed class) and outburst intensity. The slowest novae (longest decline times) occupy 
the domain of weakest outbursts and vice versa. 

In Table 4 the range of variation of some of the main nova properties is displayed, 
according to the new extended grid of results (for WD-mass range of 0.65 — 1.00 Mq). 
The parameter combination for which each maximum and minimum value was obtained is 
indicated next to the value. We note that most of the maximum values result from the new, 
lower M, runs. 

In Table 5 we present observed ranges of several nova characteristics (cf. Prialnik 1995) 
together with the ranges obtained from the full grid, for comparison. The ranges resulting 
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from model calculations overlap and completely cover the observed ranges. 


3. THE CONSTRAINED PARAMETER SPACE 

As the extension of the three-dimensional parameter space for novae has resulted in 
the emergence of new features, the question arises: where do we stop, if at all, extending 
each one of the parameters’ ranges? In other words, is the parameter space limited ? We 
examine this question based on analytical considerations, verihed by the results of numerical 
modeling. 


3.1. Heating versus cooling 


The characteristic timescale for cooling of a WD, TcooI, is basically a function only of the 
WD temperature (Mestel 1952). On the other hand, the accretion timescale Tacc is directly 
determined by the accretion rate and indirectly affected by the other parameters through 
their influence on the mass required to ignite hydrogen, Tacc = T^acc{MwD,Tw£i, M)/M. A 
nova outburst can take place only if Tcooi > Tacc ; otherwise the temperature cannot rise to 
the point of thermonuclear instability. A limiting surface is thus obtained by equating these 
timescales: 


Tcooi{Twd) — Tnacc{MwD,TwD, M)/M = 0. (1) 

Analytically, following Mestel (loc. cit.), if we use the Kramers opacity law for the WD 
atmosphere, the cooling timescale of a CO WD may be approximated by 


Tcooi ~ 2.5 X 10® 
while the core temperature is approximated by 

Twd ~ 4 X lO'^ 

On this approximation, then. 


V L/Le ) 


yr, 


( l/l^ 

VM/Me 


2/7 


K. 


Tcooi{Twd) ~ 8 X 10' 


/10^K\ 


2.5 


yr. 


( 2 ) 


( 3 ) 


( 4 ) 


\Twd J 

An approximation to the accreted mass can be derived from the requirement that the pressure 
at the base of the envelope 


Ph = 


GMwDTTlacc 
At jd) 


( 5 ) 
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exceed a critical prescribed value, say, Peru ~ 10^® dyn cm“^, in order to obtain a TNR 
(Fujimoto 1982). Fig. 3, based on grid results, presents yet another numerical validation of 
Peru having values of the order of 10^® — cm~‘^. We note that at low WD temperatures 
the range of Peru values is wider, the critical pressure changing with both M and 
Finally, using the Nauenberg approximation (Nauenberg 1972) to the Hamada-Salpeter 
WD mass-radius relation (Hamada & Salpeter 1961), 


Rwd 

Rq 


^ 1.12 X 10“^ 


/ Mwd 

V ^Ch 


- 2/3 


MwdV'^ 
Mch ) 


( 6 ) 


we obtain rriacc = T^acciMwn)- Substitution of the analytic relations into Eq. 1 yields the 
desired surface. 


A more accurate result is obtained numerically. First, we calculate cooling curves Tivuit) 
for all the WD masses of the grid (indeed, they are almost independent of Mwd)- We dehne 
the cooling timescale at any point by TcooI = {d\nT/dt)~^ and calculate it for the Twd values 
of the grid. Then, for each pair (Twd,M), we solve macc{MwD,TwD, M) = Mtcooi(Twd) to 
obtain Mwd, by interpolation on the grid results. We thus obtain the points that satisfy 
Eq. 1, dehning a surface within the parameter space. It is the bottom surface shown in 
Fig. 4. The surface exhibits a monotonic decline in Mwd for each Twd, from lower to higher 
M, but also some decline of Mwd at the same accretion rate from higher to lower Twd- As 
mentioned earlier, the cooling timescales increase with decreasing temperatures. Because 
longer TcooI (relative to Tacc) is what we require for a nova outburst to develop, the WD mass 
restriction is weaker for lower temperatures. In principle, if rriacc were mainly a function of 
Mwd, the condition for a nova outburst to occur would be eventually satisfied for any WD 
and accretion rate, because TcooI would increase, while Tacc would remain fixed. However, 
calculations show that as Twd decreases, rriacc becomes ever more strongly dependent on 
Twd, increasing steeply as the WD cools. Thus, below a certain value of Twd, the relation 
TcooI < Tacc IS maintained while both timescales increase (cf. Schwartzman et ah (1994)). 
In other words, when mass accumulates slowly as the WD cools down, the accreted mass 
remains below the critical value while both increase with time. 


The cooling time constraint is clearly illustrated in Fig. 5, where examples are shown of 
the evolution of two characteristics: the maximum temperature within the burning envelope 
Tmax and the total luminosity L, for two parameter combinations. A range of accretion rates 
is displayed, following the evolution over a period of one full cycle, from the beginning of 
accretion through outburst and decline. The lowest accretion rate for which a ’’true” mass- 
ejecting nova is obtained is M = 5 x 10“^^ Mq yr“^ (and even then, not for all parameter 
combinations, as seen in the results grid). We see in Fig. 5 that for a slightly lower accretion 
rate value 2.5 x 10“^^ Mq yr“^, the WD cools down faster than it is able to accrete sufficient 



mass to reach TNR conditions. This is clearly seen in the Tmax{t) plot: the curve gradually 
declines, without having reached an outburst. Equally, the luminosity curve remains flat for 
the lowest M. The curves for the higher accretion rates exhibit outburst events where T^axit) 
exceeds 10® K and the bolometric luminosity increases above 4 — 5 x 10^ Lq. Note also the 
obvious but nicely demonstrated fact that as we move from lower to higher accretion rates, 
the outburst occurs earlier on the timeline, hence the recurrence period becomes shorter. 


3.2. Nuclear versus gravitational energy 


The source of energy for classical nova outbursts is nuclear energy released during the 
TNR, by burning a fraction / of the hydrogen content of the accreted mass. Thus 

Enuc f ^^^accQ 1 ('f) 

where: X ^ 0.7 is the hydrogen mass fraction in the outer layers of the nova companion 
star, and Q ~ 6 x 10^® erg g“^ is the energy released per gram of burnt hydrogen. 


The greatest part of the energy released at outburst is used in lifting the ejected shell 
from the gravitational potential well of the white dwarf. This gravitational energy may be 
roughly approximated by 


Egrav 


GM\^£)Tn 




R 


WD 


( 8 ) 


Obviously, a mass ejecting outburst can take place only if E, 
surface may be defined by requiring Enuc = Egrav, 


nuc > Egrav Therefore, a limiting 


macc{MwD, M,Twd) G Mwd _ ^ 
rnej^My/M,T\/ yd) f^Q Rwd 

As a rough analytical estimate, the difference between niacc and mej, may be neglected, in 
which case Eq.(9) simply imposes an upper limit on for a given value of / (e.g., 1.22 Mq 
for / = 0.1 and 1.42 Mq for / = 0.3), a flat surface parallel to the \Ty/j:,,M] plane in the 
parameter space. This is already a more severe constraint than just the Chandrasekhar mass 
limit. Using the results of the numerical computations and taking / = 0.3 for illustration, we 
obtain a more significant constraint in the form of a slightly curved surface, the top surface 
in Fig. 4. 


3.3. Accretion versus Eddington luminosity 

For the nova progenitor to be able to accrete material during the quiescence phase, 
the accretion luminosity Lace must be lower than the Eddington critical luminosity LEdd- 
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Otherwise, radiation pressure would push away and dissipate the accreted material. In fact, 
the total (net) luminosity of the accreting star should be lower than the Eddington limit, but 
in most cases the WD intrinsic luminosity is negligible compared with Lace- The accretion 
luminosity is given by 

Lacc = OlGMwdM / RwDi ( 10 ) 

where a ~ 0.15, taking into account accretional heating (Regev & Shara 1989), and again, 
using Eq. 6 for Rwd{Mwd)- Thus a third limiting surface, obtained by equating Lace and 
LeM) 

. Attc 

M — - Rwe{Mwe) = 0 (11) 

aKs 

(where Kg is the electron scattering opacity coefficient) is shown in Fig. 4. We note that the 
result, which in this case relies solely on analytical considerations, is independent of the WD 
temperature, although it might be indirectly affected by it to some extent through a more 
realistic opacity coefficient. It is not surprising that the allowed accretion rates decrease 
with increasing Mwd- 


3.4. Degeneracy and the WD core temperature 


According to Fig. 4, the relevant region in parameter space where nova eruptions can 
occur is tube-shaped. Our calculations assumed that the relevant range of WD core temper¬ 
atures is 10 — 50 X 10® K, taking the discrete representative values of 10, 30, and 50x10® K. 
In fact, the WD temperature is restricted from above and below, as we shall now show. 


We may obtain a rough estimate for the upper limit on Two by setting the Fermi 
parameter £/ oc (InP —2.51nT) to zero at the base of the accreted layer, where the pressure, 
approximately given by Eq. 5, is of the order of 10^® dyn cm“^. Thus assuming that T^e ~ 
Th, we have 


Two ~ 4.7 X 10^ 


/ Pe/10^®dyn cm ^ 

V 



( 12 ) 


where Pe is the electron pressure. Hence if the condition of strong electron degeneracy is 
expressed as ^ 0, it results in the requirement ^ 5 x 10^ K. Therefore the tube of 
Fig. 4 has a high temperature end at a T^e of about 50 million degrees. Does it have a low 
temperature end? 


A lower limit on Two results from the following considerations. When material starts 
accumulating on the WD surface, its temperature is lower than the core temperature and 
also lower than the ignition temperature (otherwise hydrogen would ignite immediately and 
quietly, rather than explosively under degenerate conditions). As the material becomes 
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compressed, it releases gravitational energy, which is absorbed, in part, by the accreted layer, 
while in part it is conducted into the core. If the absorption of heat is sufficiently effective 
to raise the temperature of the hydrogen-rich material then, eventually, the temperature at 
the bottom of the hydrogen-rich layer will become high enough for hydrogen to ignite. The 
nuclear luminosity, low at first, will soon become the dominant energy source. Therefore, the 
restrictive condition for a nova outburst to occur is that compressional heating be sufficiently 
effective in order to raise the temperature to the ignition value of roughly 15 x 10® K required 
by the CNO cycle. The competition between compressional heating power and the rate 
of heat conduction can only be decided by solving the energy equation. Although it is 
reasonable to assume that accumulation of matter at a high rate and/or a massive WD 
will tip the scale in favor of heating, an accurate analytical result is difficult to obtain; it 
would require too many simplifying assumptions and approximations. We therefore resort 
to numerical modeling in order to obtain an estimate on the limiting temperature. 

A few additional models, beyond the grid described in Section 2 were calculated for 
WDs cooler than 10^ K; indeed very old WDs in binary systems may reach temperatures 
below 10^ K (Nelson et al. 2004). For high accretion rates (M = IO^^Mq yr“^), we obtained 
nova outbursts even for = 3 x 10® K and the results were very similar to those obtained 
for Two = 10’^ K. At the other M end, adopting M = yr“^ and Mwd = 1 Tf©, 

we did not obtain an outburst for Tvcd/ 10® K=5, 7 and 8; the accreted material cooled 
continuously. However, for Two = 9 x 10® K a regular nova outburst was obtained. Thus 
the lower limit for Twd is strongly dependent on M, as expected, and in agreement with the 
conclusions of Schwartzman et al. (1994). In the case of a very cold WD, heating is impeded 
by heat conduction into the hydrogen-depleted deep core and may be altogether suppressed. 
However, the impediment is less severe at high accretion rates, due to the rapid supply of 
accretion energy. 

We also checked the effect of Mwd'- taking Twd/^^^ K=8 (no nova for Mwd = 1 Mq), 
but increasing the WD mass to 1.25Mq, we did obtain a nova outburst. In this case we 
found the lower limit to be below 8 x 10® K but above 7 x 10® K. Thus, a lower limit to the 
WD core temperature exists for any combination of WD mass and accretion rate, but for 
high accretion rates this limit is so low that it requires more than the age of the universe for 
a cooling WD to reach it. Hence the constraint on Twd becomes signihcant only in the case 
of low accretion rates. 

Finally, since it is not clear that erupting WDs may cool down to the theoretical lower 
limit, and since between the lower limit and 10^ K the results are not highly sensitive to the 
precise value of Twd (cf. Schwartzman et al. (1994)), we chose 10^ K as the lowest Twd 
value in the present study, as in Paper I. However, accretion on very cold WDs will be further 
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investigated in a future paper, since it is this part of the parameter space that appears likely 
to shed light on the most peculiar observed novae, showing very high luminosities, relatively 
low expansion velocities and, apparently, large ejected masses. 


3.5. The full parameter space 

The three surfaces obtained so far and displayed in Fig. 4 describe a restricted, tube¬ 
shaped region within the 3D parameter space where conditions for nova outbursts are ex¬ 
pected to be satished. The additional constraints imposed by the electron degeneracy and 
mainly concerning the WD temperature determine the ends of the tube. 

Hence there appears to be a conhned volume of parameter space, where conditions 
for classical nova outbursts are satished. Plotted on top of these surfaces are all of the 
parameter combination positions that produced ’’well-behaved” mass-ejecting nova outbursts 
in our runs (including the three runs for the lowest WD mass of O.4OM0 and two runs for 
WD temperatures lower than 10^ K as analyzed in the previous section). There are six 
’’mass ejecting” parameter combinations that he on the boundaries of the restricting volume, 
denoted in Fig. 4 by black circles, while the other 75 combinations he well inside the volume, 
and are denoted by blue asterisks. As we have explained earlier, the conhning volume we 
have constructed emphasizes the existence of limitations (on each parameter range) and also 
a fairly clear shape, but one that should be understood as qualitative. Positions of mass- 
ejecting nova parameter combinations of the numerical calculations may lie on the verge or 
even slightly outside that volume, since construction of the different surfaces involves many 
constants, coefficients and approximations, influencing the exact values obtained. It would 
have been very instructive were we also able to place on this graph the positions of (parameter 
combinations of) observed novae. Unfortunately, the estimation of all three parameters for 
any observed eruption is still problematic and involves too many uncertainties. A unique 
estimation of the three basic parameters that lead to a set of observed characteristics is an 
ambitious aim yet to be achieved; hopefully it will be, with the help of studies such as the 
present one. 


4. SUMMARY AND CONCLUSIONS 

In the framework of this study we have extended the 3D parameter grid of multicycle 
nova evolution to cover combinations with low accretion rates. The main results of this study 
may be summarized as follows: 
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1. The entire range of observed nova characteristics is thoroughly covered by the com¬ 
plete grid of models. Even exceptional observed values such as outburst amplitude of over 
19 magnitudes and very high Z are covered by the new grid of results. The majority of the 
calculated maximum values for the various characteristics is obtained for the lower-M runs, 
down to the lowest accretion rate value of 5 x 10“^^ Mq yr~^. 

2. There has been a problem deriving ejected mass values as high as those deduced 
observationally. Moreover, a strong case has been made that much of the ejecta from each 
nova is missed by observers (Ferland 1998). We have managed to produce a notably high 
rriej value of 6.6 x 10“^ Mq, which is almost a factor 3 greater than the maximum value 
obtained in the original grid (Paper I). This high niej was produced for a low Mwd at low M 
(combination 065.10.12). The former maximum rriej was obtained for the same 
combination but with M of 10“^° Mq yr“^. 

3. For the new runs of lower M values we obtain maximum outburst luminosities 
surpassing the Eddington luminosities (calculated for electron scattering opacities) by factors 
of up to a few tens, for the whole range of WD masses. The highest super-Eddington 
luminosities (with correspondingly highest derived outburst amplitudes) are obtained for 
fast to very-fast novae. 

4. We predict the existence of remarkably small amplitude novae with decline times 
ranging from a week to over a year. 

5. In previous studies, as well as the current one, we have seen that by extending 
the parameter space new features emerge. That raised the question whether the parameter 
space is at all limited, and if so, what are the exact limits on the three basic independent 
parameters? Using both models and analytic relations, we derived a conhning volume within 
the three dimensional parameter space, where conditions for nova outbursts are satished. 
Thus the parameter space is indeed limited. As the new grid of models now covers the entire 
parameter space (assuming WDs composed of CO — but see below), it should enable us to 
provide explanations for most (all?) novae, including the more peculiar objects among them. 
This, however, will be the subject of a later study. Of special signihcance is the existence 
of a critical accretion rate Mcrit{Mwo-,Two) below which conditions are not sufficient to 
trigger a TNR. The critical mass transfer rate is indeed very low, but seems to be relevant 
for old close interacting binaries with very low mass companions. For cases like these, the 
determination of limits is important. 

6. Although this study covers the entire range of WD masses, practically up to the 
Chandrasekhar limit, the assumed composition of the WD core is C and O (in equal mass 
fractions) for all models. At hrst glance, this might appear unrealistic; very massive WDs are 
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believed to emerge after the carbon burning stage in the evolution of a star, and they should 
be composed of oxygen, neon and magnesium (Gutierrez et ah 1996). Although the mass 
fractions of these elements and their dependence on the WD mass are still uncertain issues, 
clearly such WDs should not contain carbon. On the other hand, our grid of models produces 
results that cover the entire range of observed characteristics (ignoring the breakdown of the 
total heavy element mass fraction Z^j). This seems to indicate that they are more widely 
applicable than imposed by the assumption on the WD composition. In fact, the mechanism 
of nova eruptions revolves around thermonuclear instability under conditions of electron 
degeneracy. Electron degeneracy, in turn, is determined by /ie =< A./Z >, which assumes 
the same value for helium and carbon burning products. In addition, the CNO cycle is 
not sensitive to the initial CNO breakdown. Consequently, we should not expect signihcant 
differences in the outburst characteristics of WDs that differ only in composition, except for 
the breakdown of Z^j. In order to test this expectation we have repeated the calculations 
for two illustrative models after replacing the carbon in the WD core by neon. The results 
are summarized in Table 6 and indeed they confirm our prediction. Thus it is not surprising 
that the entire range of nova characteristics is reproduced by our CO models and the results 
should be applicable to novae in general. 

7. An intriguing issue related to the investigation of novae is the question of the ultimate 
fate of the WD. Will it be losing or gaining mass after multiple successive cycles? Will it be 
able to reach Mchi thus acting as a possible SNIa precursor, or not? The ratio mej/niacc is 
shown in Fig. 6 for all parameter combinations. We note that it falls below unity only in a 
small region of the parameter space, thus strongly reducing the possibility of SNIa to result 
from accreting WDs. Nevertheless, this region does lead continuously from low mass to the 
Chandrasekhar limit provided the accretion rate remains very high all along. Therefore, it is 
possible, at least in principle, for a WD to grow by accretion up to Mch- In order to further 
investigate this problem we need to consider the dynamic evolution of the binary system 
that determines the evolution of mass transfer rate. Given the accretion rate as a function of 
time Macc{t), which changes with binary separation and masses, we can calculate the change 
of Mwd with time, 

Mivoit) = J [Maccit) — Mej{Macc{t), Mwiy(t))]dt, (13) 

where 

Mej{Macc{t),MwD{t)) = Ulej I Free (14) 

is obtained by interpolation on the grid. An example is shown in Fig. 7, based on an evolving 
accretion rate kindly supplied by Jarrod Hurley, for a binary system of initial masses 0.6 Mq 
(secondary) and 0.9 Mq (primary). 
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8. Finally, in a similar manner, the evolving luminosity of an accreting WD may be 
constructed as a quasi-periodic function of the variables Lace, L^ax, hMi ^rec supplied 
by the grid, all of which change with time as and Mace change. Thus, in future, an 

important use of the nova grid will be the parameterization of nova evolution for long term 
modeling of large stellar systems. 


We thank Jarrod Hurley for providing us with his results on the genesis of low-mass white 
dwarfs in close binaries (the models of Table 1 and Fig. 7). We also thank an anonymous 
referee for valuable comments and suggestions. 
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Table 1. 


Example of Binary Evolution Sequence - Creation of a Low-Mass CO WD 


Time (Myr) 

Ml 

M2 

Types 


Sep. (AU) 

Rl/RLl 

R2/RL2 


0.000 

2.809 

0.196 

MS 

MS 

110.363 

0.03 

0.01 

initial 

450.482 

2.809 

0.196 

HG 

MS 

110.363 

0.07 

0.01 

ev change 

453.303 

2.809 

0.196 

GB 

MS 

110.408 

0.19 

0.01 

ev change 

456.049 

2.808 

0.196 

GB 

MS 

49.259 

1.01 

0.03 

begin RLOF 

456.049 

0.414 

0.196 

HeMS 

MS 

0.854 

1.01 

0.03 

common env 

758.047 

0.413 

0.196 

HeGB 

MS 

0.782 

0.25 

0.92 

ev change 

786.929 

0.412 

0.196 

COWD 

MS 

0.765 

0.05 

0.94 

ev change 

1010.000 

0.412 

0.196 

COWD 

MS 

0.718 

0.05 

1.01 

begin RLOF* ** 

12000.000 

0.412 

0.040 

COWD 

MS 

0.841 

0.03 

1.06 

max time 


*MS = Main Sequence, HG = Hertzsprung Gap, GB = Giant Branch. 

**From this stage onwards CV behavior - continuous mass-loss from secondary to binary with 
periodic mass ejections. (RLOF = Roche-lobe overflow) 
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Table 2. Complete Grid Results - Characteristics of The Nova Envelope 


Param. Comb. Nova Envelope Characteristics 


Mwd 

(Me) 

Two 
(10® K) 

logM 

{Mq yr-^) 

(Mq) 


rriei 

(Mq) 

Penn 

V ■ 

J ej 

^env 

Zej 

'^8,max 

(10® °K) 


10 

-9 

A.12E- 

04 

4.55E - 04 

0.2568 

0.2579 

0.1154 

0.1204 

0.90 


10 

-10 

5.62E- 

04 

6.76E - 04 

0.2358 

0.2364 

0.1832 

0.1893 

0.98 

0.40 

10 

-11 

5.S7E- 

04 

6.96E - 04 

0.2444 

0.2322 

0.1539 

0.2075 

0.95 

0.65 

10 

-6 

8.S5E- 

06 

O.OOE + 00 

0.2830 


0.0208 


0.81 

0.65 

10 

-7 

2A5E- 

05 

O.OOE + 00 

0.3434 


0.0203 


1.10 

0.65 

10 

-8 

I.OIE- 

04 

1.03E - 04 

0.3754 

0.3865 

0.0207 

0.0215 

1.34 

0.65 

10 

-9 

1.61E- 

04 

1.63E - 04 

0.2540 

0.2650 

0.1201 

0.1350 

1.39 

0.65 

10 

-10 

2.55E- 

04 

2.76E - 04 

0.2407 

0.2498 

0.1637 

0.1786 

1.67 

0.65 

10 

-11 

2.58E- 

04 

2.37E - 04 

0.2570 

0.2538 

0.1233 

0.1847 

1.39 

0.65 

10 

-12 

3.9AE - 

04 

6.66E - 04 

0.2660 

0.1799 

0.0500 

0.4458 

1.59 

0.65 

10 

-12.3* 

5.40S - 

04 

5.04E - 04 

0.2644 

0.2723 

0.0504 

0.0733 

1.51 

0.65 

30 

-6 

8.63E- 

06 

O.OOE + 00 

0.2830 


0.0208 


0.80 

0.65 

30 

-7 

2.5AE- 

05 

O.OOE + 00 

0.3690 


0.0203 


1.09 

0.65 

30 

-8 

1.02E- 

04 

I.OIE - 04 

0.3689 

0.3801 

0.0208 

0.0214 

1.22 

0.65 

30 

-9 

l.llB- 

04 

1.21E - 04 

0.2884 

0.2956 

0.1037 

0.1102 

1.21 

0.65 

30 

-10 

9.55B- 

05 

1.21E - 04 

0.2462 

0.2489 

0.2247 

0.2370 

1.24 

0.65 

30 

-11 

5.96B- 

05 

9.46E - 05 

0.1804 

0.1817 

0.3767 

0.3910 

1.21 

0.65 

30 

-12 

4.12B- 

05 

9.61E - 05 

0.1203 

0.1216 

0.5961 

0.6085 

1.28 

0.65 

30 

-12.3 

4.44B - 

05 

1.07E - 04 

0.1149 

0.1161 

0.6170 

0.6300 

1.32 

0.65 

50 

-6 

8.94E - 

06 

O.OOE + 00 

0.2860 


0.0206 


0.80 

0.65 

50 

-7 

2.66B- 

05 

O.OOE + 00 

0.3497 


0.0203 


0.97 

0.65 

50 

-8 

1.06S- 

04 

9.88E - 05 

0.3669 

0.3774 

0.0207 

0.0212 

1.31 

0.65 

50 

-9 

7.41B- 

05 

9.1QE - 05 

0.2556 

0.2584 

0.2098 

0.2208 

1.16 

0.65 

50 

-10 

5.23E- 

05 

6.72E - 05 

0.2549 

0.2570 

0.2442 

0.2551 

1.20 

0.65 

50 

-11 

3.86E- 

05 

5.36S - 05 

0.2480 

0.2550 

0.3100 

0.3160 

1.09 

0.65 

50 

-12 

4.54B - 

05 

l.llE - 04 

0.1320 

0.1347 

0.6049 

0.6171 

1.33 

0.65 

50 

-12.3 

1.16B- 

04 

2.61E - 04 

0.1288 

0.1472 

0.6262 

0.6191 

1.60 

1.00 

10 

-6 

2.05B- 

06 

O.OOE + 00 

0.3260 


0.0204 


1.05 

1.00 

10 

-7 

8.96B- 

06 

7.21E - 06 

0.3421 

0.4187 

0.0210 

0.0223 

1.41 

1.00 

10 

-8 

2.06B- 

05 

2.22E - 05 

0.3044 

0.3314 

0.0937 

0.1002 

1.59 

1.00 

10 

-9 

4.66S- 

05 

5.18E - 05 

0.2676 

0.2969 

0.1200 

0.1270 

1.86 

1.00 

10 

-10 

8.40S - 

05 

9.72E - 05 

0.2427 

0.2744 

0.1553 

0.1620 

2.10 

1.00 

10 

-11 

8.72E- 

05 

l.OOE - 04 

0.2475 

0.2730 

0.1529 

0.1917 

2.09 

1.00 

10 

-12 

9.28E- 

05 

1.62E - 04 

0.2429 

0.1901 

0.0332 

0.4770 

2.10 

1.00 

10 

-12.3 

l.OAE- 

04 

2.00E - 04 

0.1860 

0.1772 

0.0136 

0.5328 

2.31 

1.00 

30 

-6 

2.10E- 

06 

O.OOE + 00 

0.3160 


0.0204 


1.04 

1.00 

30 

-7 

8.7AE - 

06 

5.74E - 06 

0.3280 

0.3535 

0.0208 

0.0215 

1.34 

1.00 

30 

-8 

2.03E- 

05 

2.19E - 05 

0.3035 

0.3300 

0.0951 

0.1010 

1.57 

1.00 

30 

-9 

2.70E- 

05 

3.10E - 05 

0.2686 

0.2940 

0.1492 

0.1590 

1.68 

1.00 

30 

-10 

2.10E- 

05 

2.70E - 05 

0.2376 

0.2587 

0.2411 

0.2558 

1.66 

1.00 

30 

-11 

1.15E- 

05 

1.99E - 05 

0.1752 

0.1884 

0.4435 

0.4656 

1.60 

1.00 

30 

-12 

9.72E- 

06 

2.17E - 05 

0.1407 

0.1665 

0.5910 

0.5971 

1.66 

1.00 

30 

-12.3 

3.21E- 

05 

5.53E - 05 

0.1607 

0.2064 

0.5588 

0.5431 

2.12 

1.00 

50 

-6 

2.15E- 

06 

O.OOE + 00 

0.3220 


0.0203 


1.03 













Table 2—Continued 


Param. Comb. Nova Envelope Characteristics 


Mwd 

(Me) 

Two 
(10® K) 

logM 

{Mq yr-l) 

fTT'acc 

(Me) 


rriei 

(Mq) 

Peni; 

Y ■ 

1 ej 

^env 

Zej 

'^8,max 

(10® °K) 


50 

-7 

8.30E - 

06 

5.05E - 06 

0.3239 

0.3470 

0.0207 

0.0210 

1.39 


50 

-8 

2.27E- 

05 

2.27E - 05 

0.3625 

0.3891 

0.0227 

0.0272 

1.58 


50 

-9 

1.62E- 

05 

1.94E - 05 

0.2718 

0.2934 

0.1836 

0.1957 

1.53 


50 

-10 

1.09E- 

05 

1.42E - 05 

0.2622 

0.2809 

0.2486 

0.2636 

1.46 


50 

-11 

7.9AE - 

06 

1.19E - 05 

0.2220 

0.2450 

0.3710 

0.3790 

1.43 


50 

-12 

2.07E- 

05 

3.27E - 05 

0.2328 

0.2756 

0.4323 

0.4210 

1.78 


50 

-12.3 

2.53E- 

05 

4.69E - 05 

0.2274 

0.2748 

0.4947 

0.4764 

1.96 

1.25 

10 

-6 

4.14B- 

07 

O.OOE + 00 

0.3184 


0.0204 


1.27 

1.25 

10 

-7 

1.92B- 

06 

1.62E - 06 

0.3194 

0.3744 

0.0217 

0.0263 

1.62 

1.25 

10 

-8 

3.67E - 

06 

3.91E - 06 

0.2983 

0.3637 

0.0834 

0.0876 

1.90 

1.25 

10 

-9 

9.27E- 

06 

1.06E - 05 

0.2562 

0.3336 

0.1417 

0.1411 

2.33 

1.25 

10 

-10 

1.91E- 

05 

2.18E - 05 

0.2473 

0.3221 

0.1413 

0.1562 

2.70 

1.25 

10 

-11 

2.97E- 

05 

3.61E - 05 

0.2372 

0.3155 

0.1933 

0.2092 

3.01 

1.25 

10 

-12 

3.22E- 

05 

3.33E - 05 

0.1889 

0.2485 

0.3696 

0.4371 

3.77 

1.25 

10 

-12.3 









1.25 

30 

-6 

3.S2E- 

07 

O.OOE + 00 

0.3299 


0.0203 


1.23 

1.25 

30 

-7 

1.96E- 

06 

1.86E - 06 

0.3837 

0.4380 

0.0224 

0.0273 

1.67 

1.25 

30 

-8 

3.SAE - 

06 

4.16E - 06 

0.3027 

0.3671 

0.0996 

0.1037 

1.92 

1.25 

30 

-9 

5.22E- 

06 

5.86E - 06 

0.2730 

0.3400 

0.1305 

0.1359 

2.05 

1.25 

30 

-10 

4.35E- 

06 

5.55E - 06 

0.2423 

0.3060 

0.2380 

0.2480 

2.05 

1.25 

30 

-11 

2.3AE - 

06 

4.26E - 06 

0.1676 

0.2202 

0.4733 

0.4931 

1.99 

1.25 

30 

-12 

2.10E- 

06 

4.25E - 06 

0.1421 

0.2124 

0.5747 

0.5725 

2.01 

1.25 

30 

-12.3 

5.79E- 

06 

9.91E - 06 

0.1858 

0.2708 

0.5548 

0.5401 

2.84 

1.25 

50 

-6 

4.16E- 

07 

O.OOE + 00 

0.3148 


0.0205 


1.26 

1.25 

50 

-7 

1.96E- 

06 

1.78E - 06 

0.3374 

0.4221 

0.0218 

0.0269 

1.67 

1.25 

50 

-8 

3.69E- 

06 

4.01E - 06 

0.3026 

0.3661 

0.1019 

0.1069 

1.89 

1.25 

50 

-9 

3.18E- 

06 

3.58E - 06 

0.2803 

0.3415 

0.1498 

0.1583 

1.86 

1.25 

50 

-10 

2.14E- 

06 

2.78E - 06 

0.2584 

0.3133 

0.2523 

0.2663 

1.78 

1.25 

50 

-11 

1.62E- 

06 

2.51E - 06 

0.2210 

0.2820 

0.3930 

0.3970 

1.76 

1.25 

50 

-12 

4.43E - 

06 

5.95E - 06 

0.3387 

0.4160 

0.3356 

0.3264 

2.12 

1.25 

50 

-12.3 

5.57E- 

06 

l.OOE - 05 

0.4400 

0.5262 

0.3346 

0.3169 

2.46 

1.40 

10 

-6 

1.81E- 

08 

O.OOE + 00 

0.3372 


0.0213 


1.64 

1.40 

10 

-7 

7.71E- 

08 

5.31E - 08 

0.3087 

0.4603 

0.0218 

0.0272 

2.09 

1.40 

10 

-8 

1.64E- 

07 

1.83E - 07 

0.2977 

0.5004 

0.1277 

0.1277 

2.52 

1.40 

10 

-9 

4.12E- 

07 

4.74E - 07 

0.2622 

0.4732 

0.1536 

0.1521 

3.07 

1.40 

10 

-10 

5.90E- 

07 

6.90E - 07 

0.2543 

0.4068 

0.1591 

0.2250 

3.62 

1.40 

10 

-11 

2.59E- 

06 

3.03E - 06 

0.2332 

0.4598 

0.2077 

0.1967 

4.70 

1.40 

10 

-12 









1.40 

10 

-12.3 









1.40 

30 

-6 

1.78E- 

08 

O.OOE + 00 

0.3466 


0.0208 


1.62 

1.40 

30 

-7 

7.94E - 

08 

5.54E - 08 

0.3073 

0.4574 

0.0217 

0.0271 

2.10 

1.40 

30 

-8 

2.02E- 

07 

2.02E - 07 

0.3600 

0.5410 

0.0242 

0.0339 

2.51 

1.40 

30 

-9 

2.64E - 

07 

3.02E - 07 

0.2817 

0.4891 

0.1488 

0.1497 

2.81 

1.40 

30 

-10 

2.11E- 

07 

2.68E - 07 

0.2509 

0.4643 

0.2338 

0.2325 

2.76 

















Table 2 


Continued 


Param. Comb. Nova Envelope Characteristics 


Mwd 

(Me) 

Two 
(10® K) 

logM 

{Mq yr-^) 

^acc 

(Me) 

rriei 

(Me) 

"^env 

Y ■ 

1 ej 

^env 

Zej 

'^8,max 

(10® °K) 

1.40 

30 

-11 

1.28S-07 

2.09E - 07 

0.1968 

0.3986 

0.4235 

0.4325 

2.65 

1.40 

30 

-12 

9.30S - 08 

1.94E - 07 

0.1189 

0.3213 

0.6100 

0.6179 

2.66 

1.40 

30 

-12.3 








1.40 

50 

-6 

1.80B-08 

O.OOE + 00 

0.3354 


0.0214 


1.62 

1.40 

50 

-7 

8.09S - 08 

5.69E - 08 

0.3089 

0.4772 

0.0218 

0.0277 

2.07 

1.40 

50 

-8 

2.02B-07 

2.02E - 07 

0.3600 

0.5385 

0.0241 

0.0336 

2.48 

1.40 

50 

-9 

1.90S-07 

2.18E - 07 

0.2927 

0.4952 

0.1535 

0.1541 

2.60 

1.40 

50 

-10 

1.21B-07 

1.64E - 07 

0.2608 

0.4678 

0.2848 

0.2838 

2.48 

1.40 

50 

-11 

6.83B - 08 

1.25E - 07 

0.2220 

0.4230 

0.4750 

0.4730 

2.40 

1.40 

50 

-12 

2.44B - 07 

3.15E - 07 

0.3772 

0.6059 

0.3348 

0.3202 

4.20 

1.40 

50 

-12.3 









logM of —12.3 stands for M = 5 x 10 
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Table 3. Complete Grid Results - Characteristics of The Outburst 


Param. Comb. Outburst Characteristics 


M^yd 

Two 

logM 

'^max 

(km s~^) 

'^avg 

(km 

^4,max 

(104 Lq) 

A 

ts,bol 

{days) 

^m — l 
(days) 

Prec 

{yr) 

BSi 

10 

-9 

234 

149 

7.03 

15.6 

l.lOE + 04 

4.96S + 01 

4.12E + 05 


10 

-10 

203 

97 

11.90 

18.2 

8.06B + 03 

6.19E + 01 

5.62E + 06 


10 

-11 

819 

462 

23.88 

21.0 

5.50S + 04 

3.87E + 01 

5.87E + 07 


10 

-6 

0 

0 

1.18 

4.2 

2.65B + 04 

O.OOE + 00 

8.35E + 00 

0.65 

10 

-7 

0 

0 

1.48 

6.9 

5.67B + 04 

O.OOE + 00 

2.45E + 02 

0.65 

10 

-8 

156 

122 

1.52 

9.4 

2.93B + 04 

1.17E + 03 

l.OlE + 04 

0.65 

10 

-9 

2590 

2150 

4.76 

13.2 

3.83E + 04 

2.64E + 02 

1.61E + 05 

0.65 

10 

-10 

4210 

2650 

13.70 

16.8 

3.23B + 04 

1.17S + 02 

2.55S + 06 

mM 

10 

-11 

1300 

623 

5.98 

18.0 

1.53E + 04 

2.76E + 01 

2.58E + 07 


10 

-12 

682 

216 

20.72 

20.7 

8.69B + 03 

l.lOS + 02 

3.94E + 08 


10 

-12.3 

543 

230 

17.34 

20.9 

8.22B + 02 

6.20E + 02 

1.08E + 09 


30 

-6 

0 

0 

1.17 

4.1 

2.82B + 04 

O.OOE + 00 

8.63E + 00 

0.65 

30 

-7 

0 

0 

1.49 

6.9 

6.01B + 04 

O.OOE + 00 

2.54E + 02 

0.65 

30 

-8 

139 

125 

1.58 

9.5 

3.22B + 04 

1.22E + 03 

1.92E + 04 

0.65 

30 

-9 

210 

156 

1.65 

12.0 

1.27E + 04 

6.79E + 02 

l.llE + 05 

mM 

30 

-10 

316 

195 

7.96 

16.2 

8.38E + 03 

4.83E + 02 

9.55E + 05 


30 

-11 

471 

189 

2.16 

15.2 

4.38E + 02 

3.74E + 01 

5.96S + 06 


30 

-12 

544 

268 

2.72 

15.4 

6.08E + 02 

2.75E + 01 

4.12E + 07 


30 

-12.3 

736 

297 

2.85 

15.5 

5.41E + 02 

2.61S + 01 

8.89E + 07 


50 

-6 

0 

0 

1.14 

4.1 

2.98B + 04 

O.OOE + 00 

8.94E + 00 

0.65 

50 

-7 

0 

0 

1.42 

6.9 

6.44B + 04 

O.OOE + 00 

2.66E + 02 

0.65 

50 

-8 

159 

130 

1.60 

9.5 

4.38B + 04 

1.17E + 03 

1.06E + 04 

0.65 

50 

-9 

340 

240 

1.73 

12.0 

7.91E + 03 

i.39E + 02 

7.41E + 04 


50 

-10 

369 

208 

1.91 

13.6 

7.51E + 03 

3.59E + 02 

5.22E + 05 


50 

-11 

416 

179 

2.17 

13.7 

6.54E + 03 

3.56E + 02 

3.86S + 06 


50 

-12 

342 

176 

2.87 

16.0 

5.52E + 03 

4.00S + 01 

4.54E + 07 


50 

-12.3 

552 

487 

25.32 

16.5 

1.56E + 03 

1.69E + 01 

2.33E + 08 

HiJ« 

10 

-6 

0 

0 

2.96 

4.2 

2.43B + 03 

O.OOE + 00 

2.05S + 00 


10 

-7 

265 

240 

3.45 

6.9 

2.72B + 03 

2.10E + 02 

8.96S + 01 


10 

-8 

351 

271 

3.26 

9.4 

1.35B + 03 

1.27E + 02 

2.06E + 03 

Mn 

10 

-9 

525 

256 

3.88 

12.0 

1.13E + 03 

9.35S + 01 

4.66E + 04 

mn 

10 

-10 

1920 

1180 

11.30 

15.7 

1.28B + 02 

3.36S + 01 

8.40E + 05 

Mn 

10 

-11 

1250 

1063 

6.08 

17.0 

2.96E + 03 

3.54E + 01 

8.72E + 06 


10 

-12 

1230 

292 

9.80 

19.1 

4.19E + 03 

5.53E + 01 

9.28E + 07 

B9 

10 

-12.3 

1070 

295 

7.46 

18.7 

7.14E + 03 

5.70E + 01 

2.08S + 08 


30 

-6 

0 

0 

2.96 

4.2 

2.53B + 03 

O.OOE + 00 

2.10E + 00 

Mn 

30 

-7 

267 

237 

3.26 

6.9 

3.35B + 03 

2.17E + 02 

8.74E + 01 

Mn 

30 

-8 

355 

274 

3.28 

9.4 

1.32B + 03 

1.57E + 02 

2.03S + 03 

Mn 

30 

-9 

475 

324 

3.61 

12.0 

1.02E + 03 

1.41E + 02 

2.70E + 04 

mn 

30 

-10 

512 

344 

3.55 

14.4 

7.08E + 02 

1.18S + 02 

2.10E + 05 

Qn 

30 

-11 

700 

446 

4.17 

15.3 

3.93E + 02 

1.03S + 01 

1.15E + 06 


30 

-12 

608 

414 

9.14 

16.5 

5.28E + 02 

1.06E + 01 

9.72E + 06 

B9 

30 

-12.3 

3520 

3490 

67.64 

18.5 

2.34B + 02 

3.82E + 00 

6.42E + 07 


50 

-6 

0 

0 

2.97 

4.3 

2.60B + 03 

OWE + 00 

2.15E + 00 
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Table 3—Continued 


Param. Comb. Outburst Characteristics 


Mwd 

Two 

logM 

'^max 

(km s~^) 

'^avg 

(km 

^4,max 

(104 Lq) 

A 

h,bol 

{days) 

^m — l 
(days) 

Prec 

{yr) 


50 

-7 

300 

239 

3.21 

6.8 

3.56B + 03 

2.53E + 02 

8.30S + 01 

Hnn 

50 

-8 

360 

272 

3.60 

9.5 

2.51B + 03 

2A0E + 02 

2.27E + 03 


50 

-9 

400 

328 

3.57 

12.0 

9.50E + 02 

1A3E + 02 

1.62E + 04 


50 

-10 

384 

302 

3.71 

14.0 

8.48E + 02 

9.14E + 01 

1.09E + 05 


50 

-11 

680 

471 

4.89 

14.3 

5.28E + 02 

1.22E + 02 

7.94E + 05 

Rn 

50 

-12 

736 

597 

4.81 

14.7 

4.65E + 02 

1.16S + 01 

2.07E + 07 


50 

-12.3 

3170 

3083 

89.91 

18.4 

2.30B + 02 

7.27>E + 00 

5.05E + 07 

1.25 

10 

-6 

0 

0 

4.57 

4.0 

3.11B + 02 

O.OOE + 00 

4.14E - 01 

1.25 

10 

-7 

414 

346 

4.84 

6.5 

2.59E + 02 

6.51E + 01 

1.92E + 01 

1.25 

10 

-8 

1110 

568 

6.38 

9.3 

1.31B + 02 

2.02E + 01 

3.67E + 02 

1.25 

10 

-9 

1480 

413 

6.67 

11.9 

9.78E + 01 

9.56E + 00 

9.27E + 03 

1.25 

10 

-10 

2230 

1940 

7.14 

14.5 

1.38E + 02 

2A3E + 01 

1.91E + 05 

1.25 

10 

-11 

4240 

3857 

69.37 

19.1 

4.32B + 02 

9.31E + 00 

2.97E + 06 

1.25 

10 

-12 

9680 

8267 

184.16 

21.6 

4.44E + 02 

5.86S + 01 

3.22E + 07 

1.25 

10 

-12.3 








1.25 

30 

-6 

0 

0 

4.52 

4.0 

2.46B + 02 

OWE + 00 

3.82E - 01 

1.25 

30 

-7 

474 

373 

4.77 

6.5 

2.36B + 02 

6.01E + 01 

1.96S + 01 

1.25 

30 

-8 

734 

531 

4.78 

9.0 

1.62B + 02 

3.22E + 01 

3.84E + 02 

1.25 

30 

-9 

1030 

589 

5.14 

11.6 

1.25E + 02 

3.14E + 01 

5.22E + 03 

1.25 

30 

-10 

1230 

678 

6.04 

14.3 

1.13E + 02 

2.20S + 01 

4.35E + 04 

1.25 

30 

-11 

2110 

1264 

6.01 

15.4 

5.86E + 01 

2.05S + 00 

2.34E + 05 

1.25 

30 

-12 

2710 

1696 

6.60 

15.8 

4.61E + 01 

2.00E + 00 

2.10S + 06 

1.25 

30 

-12.3 

4260 

3668 

89.08 

18.8 

2.79B + 00 

2.20E + 00 

1.16E + 07 

1.25 

50 

-6 

0 

0 

4.53 

4.0 

3.15S + 02 

O.OOE + 00 

4.10E - 01 

1.25 

50 

-7 

436 

355 

4.75 

6.5 

2.55B + 02 

4.75E + 01 

1.96S + 01 

1.25 

50 

-8 

807 

533 

4.61 

9.0 

1.58B + 02 

4.30E + 01 

3.69E + 02 

1.25 

50 

-9 

689 

485 

4.72 

11.5 

1.54E + 02 

2A7E + 01 

3.18E + 03 

1.25 

50 

-10 

921 

542 

4.97 

14.1 

l.OlE + 02 

3ME + 01 

2.14E + 04 

1.25 

50 

-11 

889 

658 

5.45 

14.3 

7.77E + 01 

2.59S + 01 

1.62E + 05 

1.25 

50 

-12 

1350 

751 

6.30 

14.6 

6.91E + 01 

2.70E + 00 

4.43E + 06 

1.25 

50 

-12.3 

3130 

1675 

12.25 

15.9 

4.45B + 01 

2.08E + 00 

l.llE + 07 

1.40 

10 

-6 

0 

0 

5.92 

3.4 

1.08B + 01 

O.OOE + 00 

1.81E - 02 

1.40 

10 

-7 

1410 

681 

6.03 

5.9 

1.19B + 01 

4.13E + 00 

7.71E - 01 

1.40 

10 

-8 

3060 

1160 

9.27 

8.8 

5.02B + 00 

1.57E + 00 

1.64E + 01 

1.40 

10 

-9 

2850 

1020 

9.23 

11.3 

4.31E + 00 

1.39E + 00 

4.12E + 02 

1.40 

10 

-10 

5270 

2020 

10.10 

13.9 

5.45B + 00 

6.78E - 01 

5.90E + 03 

1.40 

10 

-11 

4660 

2817 

37.54 

17.5 

4.71B + 01 

2.95E + 00 

2.59E + 05 

1.40 

10 

-12 








1.40 

10 

-12.3 








1.40 

30 

-6 

0 

0 

5.86 

3.3 

l.llB + 01 

O.OOE + 00 

1.78E - 02 

1.40 

30 

-7 

1160 

723 

5.99 

5.9 

1.19B + 01 

3.64E + 00 

7.94E - 01 

1.40 

30 

-8 

1760 

840 

6.14 

8.4 

9.35S + 00 

4.63E + 00 

2.02E + 01 

1.40 

30 

-9 

3760 

1800 

6.98 

11.0 

6.17E + 00 

1.29E + 00 

2.G4E + 02 

1.40 

30 

-10 

4490 

2120 

7.67 

13.6 

4.20E + 00 

1.50E + 00 

2.11E + 03 
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Table 3—Continued 


Param. Comb. Outburst Characteristics 


Mwd 

Two 

logM 

Rmax 
{km s~^) 

'^avg 

{km 

^4,max 

(104 Lq) 

A 

ts,bol 

{days) 

^m — l 
{days) 

Prec 

{yr) 

1.40 

30 

-11 


2738 

8.41 

15.3 

3.57E + 00 

1.87E - 01 

1.28E + 04 

1.40 

30 

-12 

3310 

2312 

9.52 

16.1 

3.05E + 00 

1.91E - 01 

9.30E + 04 

1.40 

30 

- 12.3 








1.40 

50 

-6 

0 

0 

5.85 

3.3 

l.llB + 01 

O.OOE + 00 

1.80E - 02 

1.40 

50 

-7 

1360 

812 

6.01 

5.8 

1.25B + 01 

4.24E + 00 

8.09E - 01 

1.40 

50 

-8 

1540 

833 

6.46 

8.4 

9.49B + 00 

4.80E + 00 

2.02E + 01 

1.40 

50 

-9 

3740 

1580 

7.18 

11.1 

6.03E + 00 

1.76E + 00 

1.90E + 02 

1.40 

50 

-10 

4350 

1870 

7.67 

13.6 

4.71E + 00 

1.17S + 00 

1.21E + 03 

1.40 

50 

-11 

3630 

1900 

7.74 

14.8 

3.45E + 00 

1.67E + 00 

6.83E + 03 

1.40 

50 

-12 

3550 

2570 

9.31 

14.9 

3.81E + 00 

3.13E - 01 

2.44E + 05 

1.40 

50 

-12.3 









Table 4. Grid Results - Maximum and Minimum Values of Nova Characteristics (for 

0.65 - I. 4 OM 0 ) 


Max Param. Comb. Min Param. Comb. 


Characteristic 

Value 

Mwd 

Twd 

logM 

Value 

Mwd 

Twd 

logM 

PPacc 

5A0E - 04 

0.65 

10 

-12.3 

6.83E - 08 

1.40 

50 

-11 

rriej 

6.6QE — 04 

0.65 

10 

-12 

5.31E-08 

1.40 

10 

-7 

Zej 

0.63 

0.65 

30 

-12.3 

0.02 

1.00 

50 

-7 

Yej 

0.61 

1.40 

50 

-12 

0.12 

0.65 

30 

-12.3 

Xey 

0.65 

0.65 

10 

-12.3 

0.06 

1.40 

30 

-12 

T 

8,max 

4.7 

1.40 

10 

-11 

1.1 

0.65 

50 

-11 

Pmax 

5270 

1.40 

10 

-10 

139 

0.65 

30 

-8 

"^avg 

3860 

1.25 

10 

-11 

122 

0.65 

10 

-8 

L^^bol^max 

90 

1.00 

50 

-12.3 

1.5 

0.65 

10 

-8 

^bol,max 

-10.2 

1.00 

50 

-12.3 

-5.7 

0.65 

10 

-8 

A 

20.9 

0.65 

10 

-12.3 

5.8 

1.40 

50 

-7 

tsfiol 

120 yr 

0.65 

50 

-8 

3.05 d 

1.40 

30 

-12 

tm—l 

67 yr 

0.65 

10 

-11 

0.68 d 

1.40 

10 

-10 

P 

^ rec 

1.08i7 + 09 yr 

0.65 

10 

-12.3 

281 d 

1.40 

10 

-7 
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Table 5. Observed vs. Calculated Nova Characteristics Ranges 


Characteristic 

Observed Ranges 

Observed Exceptions 

Calculated Ranges 

^max 

—6 to —9 

-10 (V1500 Cyg) 

-5.7 to -10.2 

A 

7-16 

19.3 (V1500 Cyg) 

5.8-20.9 

0 

4 — 300 days 

SymN:more 

0.76 d-Q7 yr 

rUej 

1 - 30(E - OSMq) 

RNdess 

5.3E - 08 - 6.6E - 04 

Zej 

0.04 - 0.41 

0.86 (V1370 Aql) 

0.02-0.63 

yej 

0.21 - 0.48 

0.1 (V1370 Aql) 

0.12-0.60 

V * 

^exp 

350 - 2500 km 

SymN:~ 100 km 

122 - 3860 


^average expansion velocities 


Table 6. Ejecta Composition - CO vs. ONe Models 



1.25, 10, -11 

CO ONe 

1.40, 10, -9 

CO ONe 

rriej 

Uej 

Zej 

3.61E-05 

0.3155 

0.2092 

3.67E - 05 

0.3123 

0.2179 

4.74E - 07 

0.4732 

0.1521 

2.89E-07 

0.5129 

0.1508 

* (^12 

3.353E - 02 

1.970E - 02 

2.430E - 02 

1.143E - 02 

^13 

3.136E-02 

2.089E - 02 

2.485E - 02 

8.148E - 03 

Ari4 

6.545E - 02 

3.086E - 02 

9.392E - 02 

4.706E - 02 


6.660E - 02 

3.321E - 02 

4.191E - 03 

1.147E - 03 

^16 

1.165E-02 

9.241E - 05 

2.092E - 04 

1.402E - 04 


6.157E - 04 

7.729E - 05 

1.342E - 05 

8.500E - 06 

Ne 

3.384E - 07 

6.579E - 02 

1.372E - 07 

1.358E - 03 

Na 

3.535E - 09 

6.541E - 04 

1.047E - 10 

1.141E - 06 

Mg 

2.861E-08 

9.611E-03 

5.238E - 09 

3.640E - 05 

A/26 

2.098E - 08 

7.929E - 03 

6.355E - 09 

4.258E - 05 

A/27 

1.027E - 08 

3.137E - 03 

1.251E - 09 

7.662E - 06 

> Si 

1.046E - 07 

2.600E - 02 

6.087E - 07 8.147E - 02 


Heavy elements Zf.j breakdown 
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Fig. 1.— Time of decline t^^vis (order of tm-i) versus the outburst amplitude A for both 
observed NA + NB galactic novae (taken from Duerbeck 1987) - circles, and our model 
results: crosses - M ^ 10“^°, asterisks - M ^ 10“^^. 
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Fig. 2.— 3D plots of four major properties out of the complete result grid. Each plot 
represents the property as function of [M, Mwd] for a representative WD temperature (stated 
in units of 10® K). In panel a, the white surface represents the critical Eddington luminosity, 
emphasizing the domains where the luminosity either surpasses or is lower than LEdd- 








































26 



Fig. 3.— Values of the pressure at the base of the envelope prior to the TNR, grouped by 
the three main values as function of the accretion rate; calculated from grid results 
(according to eq. 5). 
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Fig. 4.— The combined restricting surfaces, confining a volume within the 3D parameter 
space where conditions for nova outbursts are satished. The bottom surface corresponds to 
the heating versus cooling criterion (section 3.1); the top WD mass-limiting surface corre¬ 
sponds to the nuclear versus gravitational energy considerations (section 3.2); the left mass 
transfer rate-limiting surface relates to the accretion versus Eddington luminosity (section 
3.3). Positions of the grid parameter combinations that successfully produced nova eruptions 
are plotted within the conhned volume. 
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Fig. 5. — Evolution in time of two characteristics for different accretion rates: the max¬ 
imum temperature (within the burning shells) T^ax and the luminosity L, for M\y£, = 
I.OOMq, Ti 4 /£) = 30 X (left), and = 1.25Mq,Tiyd = 50 x 10®iF (right). 
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Fig. 6.— Color maps displaying shades that correspond to values of the ratio rriej/macc in 
the [Mwd, LogM] plane. A value of 1 corresponds to an unevolving WD mass, below - an 
increase in Mwd, and above - a decrease in time of Mwd- The three panels correspond to 
the three values - 10, 30 and 50 (lO^iF), top to bottom. 
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Fig. 7.— An example calculation of M\Y£,{t) based on an evolving M for initial masses of 
0.6 and 0.9 Mq of the secondary and primary respectively. 







